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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE NO • 1563 


EFFECT OF PARTIAL WING- LIFT IN 
SEAPLANE LANDING IMPACT 
By Stanley IJ. Benscoter 


SUMMARY 

A solution, is presented for the motion of a prismatic float of 
infinite length that is dropped vertically into the water at zero trim. 

A lift force, simulating seaplane wing lift, is assumed to act on the 
float and to remain constant during the imp act period- The float is 
assumed to have a uniform mass per unit of length and a uniform wing lift 
per unit of length. The solution is determined for values of the lift 
ranging from zero to full wing lift. Time variations of the acceleration 
have "been computed for a specific numerical example • The value of the 
mass ratio (ratio of apparent water mass to float mass) at the instant of 

maximum acceleration is determined in general form. The variation _of 

maximum acceleration with sinking speed is i 3-lustrated for various 
amounts of wing lift. 


INTRODUCTION 


During a rough-^water landing of a seaplane, some stalling of the 
wings may occur during the second and subsequent impacts. The wing lift 
becomes less than the weight of the seaplane. Solutions for the maximum 
acceleration in smooth water landings are available for both two— and 
three-dimensional cases of fluid flow when the wing lift equals the 
weight of the seaplane (references 1 and 2). These solutions may be 
converted to rough-water landings by assuming a reasonable shape for the 
wave surface. With partial wing lift the solution for the three— 
dimensional case becomes very complex mathematically. The effect of 
partial wing 3J_ft upon the maximum acceleration occurring during the 
impact is not known. However, the two-dimensional case can be readily 
solved. This solution offers valuable suggestions- for the writing of 
structural design specifications. 


SYMBOLS 


a acceleration 

c half loaded width 

dimensionless parameter 


C 
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g acceleration due to gravity 

H momentum of fluid 

lc wing lift factor (L/W) 

L wing lift per unit length 

m apparent water mass 

M float mass per unit length 

t time 

W weight of- float per unit-length 

z draft of float 

0 dead-rise angle 

e apparent-mass coefficient 

T] dimensionless parameter 

p. mass ratio ( m/M ) 

Subscripts: 

o instant of entry (t = 0) 

m instant of maximum acceleration 

A dot is sometimes used to indicate differentiation with respect to 
time. - 

The magnitude of all physical quantities is assumed to he determined 
by a consistent system of units. 

SOLUTION OF DIFFERENTIAL EQUATION 


In the following analysis, a two-dimensional solution of a float 
entering water is made] that is, a prismatic float of infinite length 
is considered to enter the water at zero trim from a vertical drop. A 
lift force is assumed to he distributed over the length of the float 
. in correspondence to the wing lift in a seaplane landing and is assumed 
to remain constant— during the impact. The float- is assumed to- have a 
uniform mass per unit of length and a uniform wing lift - per unit of 
length. A croBS section of the float and coordinate system used is 
shown in figure 1. In the analysis, the equations are written for a 
unit slice of float and of fluid. The float has a mass M per unit 
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length. The fluid is assumed to have a total momentum that can "be 
defined in terms of an apparent mass. Xn the present case^ the apparent 
mass may "be considered to he the mass of fluid contained in a had f circular 
cylinder as originally introduced hy von Karman (reference 1) and shown 
in figure 1. The float is acted upon hy a lift force ~L, its own 

weight ¥, and the water reaction. Forces and drafts are considered ~~ — 

to he positive downward. 


The apparent mass is proportional to the square of the loaded wi dth. 
Consequently , with flat— side wedge— bottom floats^ the apparent mass is 
^■lso proportional to the square of the draft* The apparent mass and 
momentum of a unit slice of the fluid may he expressed in the following 
forms: 


m = 

> 

H = mz 


( 1 ) 


In order to write an equation of motion governing the draft of the 
float , the float and the fluid may he considered to form a single system 
that is acted upon hy the external forces ¥ and L. Newton* s law on 
the rate of change of momentum may then he written for the system as 
follows : 


— (Mz + mz) = ¥ — L 
dt 


( 2 ) 


This equation could he developed hy a somewhat different argument. The 
second term ^ (mz) represents the rate of change of momentum of the 

fluid and thus may he regarded as the force exerted on the fluid hy the 
float. If this quantity were preceded hy a negative sign it would "become 
the force exerted hy the fluid on the float and would he written on the 
right-hand side of the equation. The right -h a n d side would then repre- 
sent an of the forces acting on the float. The equation would then he 
regarded as an application of Newton* s law of motion to the float itself. 


A fairly large amount of experimental evidence is available to 
Justify the use of this ideal fluid theory. Test results obtained in 
the NACA impact basin have shown good agreement with theory. The problem 
of the three-dimensional case for full wing lift has been solved and 
compared with experiment in reference 2. 
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The Initial conditions for the solution of- the differential equation 
are as follows: when t = 0, 


z = 0 
z = z D 


> 


( 3 ) 


= a 


o 


The initial Telocity and initial acceleration are assumed to he known. 

The lift may he expressed hy the product kW where k is considered 
herein to. vary from 0 to 1. Introducing this value into equation (2) and 
performing the indicated differentiation gives 


(M + m)z +-z-— = (1 — k)W 
dt 


00 


This equation is nonlinear hut may he readily converted into an inte— 
grahle form. It is convenient in developing the solution of this equa- 
tion to introduce the dimensionless mass ratio u and the initial 
acceleration a Q : 





( 5 ) 


a o = l o = (1 “ k ) g 


(6) 


Dividing through equation (h-) hy M and inserting the formulas ( 5 ) 
and (6) gives 


(1 + 


u) z Si = 

dt dt 


(T) 


This equation may he integrated to obtain the following formula for 
velocity: 


'o +- a o t 


(8) 


z 


1 + u 
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Tn this equation the appropriate value of the constant of integration 
has "been introduced. 

The equation for velocity may again he integrated to obtain the 
relationship between draft- and time* Equation (8) may be written in 
the following form: 

U - 4) dz = (z Q + a Q t) dt (9) 

After substitution of equations ( 1 ) and (5 ), equation (9) may be inte- 
grated to give the following formula: 

z ( 1+ ^) = *o t + ^ 


Tn this case the constant of integration is zero. This equation is a 
cubic equation in z or a quadratic equation in t. Hence, it is more 
convenient to solve for t in terms of the draft. Therefore, 



Tn the design of a seaplane the quantity that is required is the 
value of the acceleration. It may be noted that the original equation 
of motion is linear in the acceleration. Consequently, it may be solved 
to obtain a formula- for acceleration in terms of drafts and velocitie s . 
Equation (7) may be solved for the acceleration to obtain 


t • 

z 


i + 4 


z 

1+4 


dfci 

dt 


a ° _ dji 

1+4 1 + 4 dz 


( 12 ) 
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Prom equations (1) and ( 5 ) it is seen that — p — is related to the 
draft "by the formula 





( 13 ) 


Differentiating equation ( 13 ) gives 


2ez = Si 

dz M z 


( 1*0 


Substituting equations ( 8 ) and (l4) into equation (12) giveB the following 
formula for acceleration: 


a o _ 2 ^ ( z 0 + a p ^) 2 
1+4 2 U+4) 3 


( 15 ) 


Equation (15) may be regarded as a formula for the acceleration in terms 
of time and draft. It is convenient to -eliminate t — in order to obtain 
a formula for acceleration in terms of 1 the draft-. — 

Equation (11) may be rearranged by transposing the first term on 
the right-hand side and squaring both sides. 


( z o + V ) 2 = *0 + 2 a o z (i + f) 


(16) 


This formula may be substituted into equation ( 15 ) to obtain 

z = - g tiz n 2 + a p(3 - 6n - P 2 ) 
z ( 1 +- p) 3 3(1 + p) 3 


( 17 ) 


If various values of z are assumed, the corresponding values of- z can 
be computed, use being made of equation (13 ) to determine p. The asso- 
ciated values of t may be computed from equation (11). A graphical 
relation between acceleration and time may then be plotted. By use of 
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formula (17) for acceleration, a formula may "be readily derived for the 
hydrodynamic reaction force. The resulting force formula would "be in 
agreement with that for the case of no wing lift ( a Q = g ) given in 

reference 3. 


NUMERICAL EXAMPLE 


For plotting purposes, it is convenient to express the acceleration 
in dimensionless foryn. Dividing through equation ( 17 ) "by g gives 


z 

g 


-2p ( " z o j + (1 - k)(3 - 6p - p 2 ) 

(1 + p) 3 ' gz ' 3(1 + p) 3 


(18) 


In order to illustrate the variation of acceleration with time a 
simple example has "been computed. The following values have teen assumed 
for the float properties and the initial conditions: 


Mass of float, M, slu^s per foot 100 

Apparent mass coefficient, 6, slugs per cubic foot 20 

Initial velocity, z Q , feet per second *12 

Wing lift factor, k 1 or 2/3 


The value k = 1 


is for full wing lift; whereas, the value 



is 


a reasonable reduction factor for wing lift to be used in design practice 
A graphical representation of the solution for the two cases is shown in 
figure 2 . It may be seen that there is some reduction in the maximum 
acceleration with partial wing lift. However, the effect is Bmall enough 
to be disregarded in meat practical design cases. 


MASS RATIO AT MAXIMUM ACCELERATION 


In design practice It Is actually necessary to’ know only the value 
of maximum acceleration that occurs during the impact. This maximum 
value may be found by differentiating the formula for acceleration with 
respect to time or draft. This derivative is set equal to zero. Thus, 


~ = 0 • ( 19 ) 

dz 
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Differentiation of equation (17) gives, 'after some rearrangement of- 
terms, - - 


-“Mm- <*)- 3(5^,-i)z 0 a - 0 (20) 


The subscript m has been introduced, to indicate that the values of 
the quantities are those that ‘occur at the instant of max imum 
acceleration. 


Equation (13 ) shovs that the draft may be expressed in terms 
of- u. Thus, 



( 21 ) 


This formula holds true at all instants, including the instant of 
maximum acceleration. Hence, it may be used to el imina te the factor 
in equation (20) by direct substitution. The result may be written in 
the following form: 


(l^ Mm ) _ ^o^ _ q 

3(5m»- 1) " a 0 \|m 


( 22 ) 


The parameter C . depends upon the initial conditions and the float 
properties. Instead of solving for it is more convenient to assume 

various values of ^ and compute the corresponding values of C. Such 

values are shown in the following table and are plotted in figure 3 . 


ha 

C 

0.20 

OO 

.22 

19.94 

.24 

10.24 

.26 

6.99 

.28 

5-35 

• 30 

4.28 

.40 

2.29 

.50 

1.53 

.60 

1.12 

.80 

.63 

1.00 

• 33 

1.20 

.11 

1.32 

0 
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Equation (22) shows that for full wing lift (a Q = 0), must 

have the value 0.2. This universal value of ^ was originally given 
in reference 3* 


MAXIMUM ACCELERATION 


A formula for the relative acceleration at any instant is given "by 
equation (18). The value of as determined from equation (22) may 

he substituted into equation (l8) to determine the maximum acceleration. 
However, it is possible to develop a simplified formula for the accelera- 
tion which is applicable only at the instant of maximum, acceleration. 

Equation (20) may be rewritten by use of equation (6), in the fol- 
lowing form: 

z£ . (1 - k)(l5 - 101^ - p, 2 ) (23) 

S=% 3(51% - 1) 

Substitution of this equation into equation (18) gives, after combining 
terms, __ 


Jm _ _ 1 — k 
S “ 1 


(2h) 


In. 
speed z 


order to show the variation of acceleration with sinking speed 
Q , the parameter C (equation (22)) may be expressed as 


where 



(25) 



(26) 


For a particular value of k, the parameter C may be determined for 
various assumed values of q from equation (25). The corresponding 
values of may be determined from equation (22) and substituted into 
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equation (24-) to obtain accelerations. These accelerations may then he 
plotted against the parameter q. Such graphs, for various values of k 
are shown in figure 4-. 

At full wing lift, equation (24-) hecomes indeterminate. The 
values k = 1 and ^ = 0.2 may he substituted into" equation ( 18 ) . 

Introducing equations (21) and (2 6 ) gives the formula 


• • 

Zjn 

g 


-25 

108\jo72 




-0.518ti 


(27) 


For this case the graph in figure 4- is a straight line. 


CONCLUSIONS 


The vertical drop of a prismatic float of infinite length with 
constant partial wing lift has heen analyzed to determine the maximum 
acceleration. The float is assumed to have a uniform mass per unit of 
length and a uniform wing lift per unit~cf length. The mass ratio at 
the instant of maximum acceleration is found to depend upon a single 
dimensionless parameter involving the initial conditions and the float 
properties. The solution shows that partial wing lift may have a small 
effect upon the maximum acceleration. The effect is small enough, 
however, to he disregarded in most practical design cases. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., February 5, 1948 
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Figure l.~ Two-dimensional floa+ 
enteri ng fluid. 
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Figure Z.~ Accelerations with -full 
and partial wing lift. 
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Figur'e 3. - Mass ra+io a+ maximum accelera+ion 
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